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Abstract: We study the spacetime obtained by superimposing two equal Aichelburg-Sexl 
shock waves in D dimensions traveUng, head-on, in opposite directions. Considering the 
colhsion in a boosted frame, one shock becomes stronger than the other, and a perturbative 
framework to compute the metric in the future of the colhsion is setup. The geometry is 
given, in first order perturbation theory, as an integral solution, in terms of initial data 
on the null surface where the strong shock has support. We then extract the radiation 
emitted in the collision by using a D-dimensional generalisation of the Landau-Lifschitz 
pseudo-tensor and compute the percentage of the initial centre of mass energy e emitted as 
gravitational waves. In D = 4 we find e = 25.0%, in agreement with the result of D'Eath 
and Payne [12]. As D increases, this percentage increases monotonically, reaching 40.0% 
in Z) = 10. Our result is always within the bound obtained from apparent horizons by 
Penrose, in L> = 4, yielding 29.3%, and Eardley and Giddings [16], in D > 4, which also 
increases monotonically with dimension, reaching 41.2% vn D = 10. We also present the 
wave forms and provide a physical interpretation for the observed peaks, in terms of the 
null generators of the shocks. 
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1 Introduction 

The ongoing physics runs at the Large Hadron Collider (LHC) with 7 TeV centre of mass 
energy, are starting to set bounds on physics beyond the standard model. For the partic- 
ular case of TeV gravity models [1], microscopic black hole formation and evaporation is 
predicted for scattering experiments with partonic centre of mass energy beyond the TeV 
[2, 3]. The first set of bounds for these models was recently released by the CMS [4] and 
ATLAS [6, 7] collaborations. The current analysis for the 7 TeV data, is however extremely 
dependent on regions of parameter space where (at best) black holes with masses close to 
the unknown Planckian regime would be produced (see [8] for a discussion). Otherwise if 
strong cuts are imposed, such that the objects produced are in the semi-classical regime 
(which is the calculable one), the cross-sections become negligible at 7 TeV, and only after 
the upgrade of the beam energy to 14 TeV, planned to take place in 2013, will the scenario 
be properly tested. Any improvement in the phenomenology of these models is therefore 
quite timely. 

Two event generators, CHARYBDIS2 [9] and BLACKMAX [10] are being used at the LHC 
to look for signatures of black hole production and evaporation. The event rates they 
produce depend sensitively on various input parameters, amongst which two fundamental 
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quantities for the modeling of the black hole production phase are: the critical impact 
parameter for black hole formation in parton-parton scattering; and the energy lost into 
gravitational radiation in this process. The latter, if large, could be an important signature 
for discovery or exclusion, if dominating the distribution of missing energy and if calcu- 
lated with enough precision in the semi-classical regime. Furthermore, if the experimental 
searches continue to rely on events with little (or under-estimated) missing energy as in [4- 
6], then this information will be crucial to determine if there is enough phase space left 
with little missing energy. 

As first argued by t'Hooft [11], classical gravity, described by general relativity, should 
dominate transplanckian scattering. Following this rationale, the best estimates, so far, 
for the two mentioned quantities, come from apparent horizon computations, initiated by 
Penrose in four spacetime dimensions. He found an apparent horizon on the past light cone 
of the collision between two shock waves, representing the gravitational field of two particles 
boosted to the speed of light. Thus, Penrose concluded that no more than 29.3% of the 
centre of mass energy could be emitted in gravitational radiation in a head-on-collision. 

This analysis was generalised to D dimensions by Eardley and Giddings [16]. They 
obtained the bound on the percentage of the centre of mass energy, eradiated, emitted in a 
head on collision, which increases with dimension and is given by 



where r2„ is the volume of the unit n-sphere. This bound increases monotonically with 
dimension approaching 50% in the limit of infinite D. Using the results in [16] for non- 
head on collisions, the critical impact parameters were numerically evaluated by Nambu 
and Yoshino [17]. According to their estimate, a black hole will form if at least half of the 
effective gravitational radius of each of the colliding particles overlap. The state of the art 
results were obtained by Yoshino and Rychkov [20] , who found an apparent horizon on the 
future light cone of the collision. Their analysis coincides with the one in [16] for the energy 
emitted in head on collisions; but the critical impact parameters for black hole formation 
become larger, yielding a larger cross section for black hole formation. Events generated 
by CHARYBDIS2 rely on a selection of random points from the configuration space allowed 
by the bounds determined in [20]. 

The event horizon, however, is located outside the apparent horizon, and the latter 
will evolve and settle to the former in the future of the collision. Thus, both the critical 
impact parameter and eradiated should be exactly determined by knowledge of the metric 
in the future of the collision ^. There are in principle two methods to compute it, both of 
which have been developed in four spacetime dimensions. 

Modeling the gravitational field of the colliding partons as shock waves, one may use 
the Aichelburg-Sexl metric [15] to describe the process. Using a method first developed by 
D'Eath [18] and D'Eath and Payne [12-14], one may compute the metric to the future of 

^Another attempt to extract the radiation using test particle interactions in Minkowski space was pre- 
sented in [22, 23]. 
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the collision analytically, as a perturbative expansion. In four dimensions, second (first) 
order perturbation theory yields 16.3% (25.0%) for the energy emitted in gravitational 
waves in a head on collision. 

Modeling the gravitational field of the colliding partons as highly boosted black holes, 
one may use numerical relativity techniques to perform fully non-linear numerical simula- 
tions of the collisions. In four dimensions the high energy collision of two black holes was 
performed by Sperhake et. al. [21], obtaining that 14±3% of the energy was converted 
into gravitational radiation. This overlaps with the result obtained by D'Eath and Payne 
in second order perturbation theory [14], which indicates that such an approximation pro- 
vides a good estimate. Moreover, introducing an impact parameter [24], one observes that 
the critical impact parameter increases from about 0.8 in [20] to about 1.2, in units of the 
Schwarzschild radius of the centre of mass energy. This is an increase of about 50% which, 
if verified in the higher dimensional case as well, leads to significantly larger cross-sections, 
thus becoming easier to exclude regions in the model parameter space. 

Phenomenologically interesting TeV gravity models occur in dimension D > 6. It 
would therefore be useful to develop each of the two above methods in higher dimensions. 
Numerical relativity for higher dimensional spacetimes is being developed (see eg. [29, 32]), 
and the first ever results for black hole collisions in higher dimensions (albeit low energy 
collisions) have been produced [30, 31]. In this paper we shall perform a higher dimensional 
collision of shock waves and determine, for even D, the energy radiated by computing the 
geometry to the future of the collision to first order in perturbation theory, building up 
on the technique of [12]. The results we obtain for the radiated energy are summarised in 
the following table, where the apparent horizon bounds obtained in [16] are also given for 
comparison: 



Spacetime dimension 


4 


5 


6 


7 


8 


9 


10 


Apparent horizon bound (%) 


29.3 


33.5 


36.1 


37.9 


39.3 


40.4 


41.2 


First order perturbation theory (%) 


25.0 




33.3 




37.5 




40.0 



Like the apparent horizon bound, our results indicate that the energy radiated increases 
monotonically with D. Moreover, our result is always below the bound, as expected. If 
this reduction trend is verified by higher order perturbation theory (or other methods, such 
as numerical black hole collisions) this result has phenomenological implications. More 
energy lost into gravitational radiation means a less massive final black hole is produced. 
Thus, this result suggests that the final black hole will be more massive than previously 
estimated, making it more consistent with the semi-classical analysis used for estimating 
the potentially observable Hawking radiation. 

This paper is organised as follows. In Section 2 we discuss the spacetime geometry of 
one and two shock waves (before the collision) in both Brinkmann and Rosen coordinates. 
In Section 3 the perturbative computation is setup and the integral solution for the metric 
in the future of the collision is provided. Some details concerning the derivation of this 
solution are provided in Appendix A. In Section 4 the extraction of gravitational radiation 
is studied. A formula in D dimensions is first derived in terms of the metric perturbation 
in the future of the collision starting from the Landau-Lifschitz pseudo tensor (in D dimen- 
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sions). Then, using results from Appendix B, a workable expression is obtained to extract 
the energy carried out by the gravitational radiation (formulas (4.7) and (4.8)). These ex- 
pressions are evaluated numerically to yield the results presented in the table above. The 
wave forms are exhibited and an interpretation is given for the domain where they have 
support and peak. We also discuss the approximations used in our derivation and further 
comment on our results in Section 5. 



2 One D-dimensional shock wave in Rosen form 

The D-dimensional Tangherlini solution with mass M is [25] 

- - - (D - 2)^0^2 )'^^ ^D- 2)nD^2 ) + r dU^., , 

(2.1) 

where dQ'j^_2 and ^d~2 are the line element and volume of the unit D — 2 sphere. The 
Aichelburg-Sexl solution [15] is found by boosting this black hole and then taking simul- 
taneously the limit of infinite boost and vanishing mass, keeping the total energy /i fixed. 
We use a coordinate system (u, f,x*), where the retarded and advanced times {u,v) are 
(t — z,t + z) in terms of Minkowski coordinates, and x* are the remaining Cartesian co- 
ordinates on the plane of the shock, i = 1 . . . D — 2, such that the transverse radius is 
p = sj x^Xi. The resulting geometry for a particle moving in the +2 direction is 

ds^ = -dudv + dp^ + p^fi$7|,_3 + K^{p)5{u)dv? , (2.2) 

where n = SirGofJ-f^D-s- The function $ depends only on p and takes the form [16] 



Hp) 



-21n(p) , D = 4 

2 (2.3) 

, D > 4 



[ {D-4)pD-^ 

The coordinates used in (2.2) are of Brinkmann type [26]. Presenting the shock wave 
in this chart, geodesies and their tangent vectors appear discontinuous across the shock. 
One can, however, introduce a new coordinate system defined by [16] 

u = u , 

V = k9{u) ^ H — — — =v + k9{u) ^ H — 



( -( 1 _i_ 9{u) J / \ 

= x' + K'^VMx)e{u) ^ } P- P_\ ^ 'ip ) ^ (2.4) 



where 6 is the Heaviside step function and $ and its derivative are evaluated at p. In this 
new chart both geodesies and their tangents are continuous across the shock at u = u = 0; 
(pa are the angles on the {D — 3)-sphere and a = 1...D — 3. Using (2.4) the metric becomes 
[18-20] 

ds^ = -dudv + (l + ^^"fdf + p\l + ^'^'fdnU , (2.5) 
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CM frame 



boosted frame (-z direction) 



Figure 1 . Spacetime diagram of the collision in the centre of mass (left panel) and boosted (right 
panel) frames. The directions are orthogonal to the t and z axis. In the centre of mass frame, the 
radiation is symmetric under z — >■ — z. Only the radiation that propagates along a small solid angle 
around = tt in the centre of mass frame (red area, delimited by the polar angle tan 6 = yT?^ — 1) 
is still propagating in the —z direction in the boosted frame. All the radiation in this solid angle is 
redshifted in the boosted frame. 



which we dub the Rosen form [27] of the shock wave. The geometry for an identical shock 
wave traveling in the —z direction is obtained by changing z — )• —z in (2.2) or equivalently 
exchange w f-)- u in (2.5). Following [12] it is simple to see that in a boosted frame (moving 
with respect to the (u, v) chart with velocity /3 in the —z direction), the oppositely directed 
shock waves keep their form, but with new energy parameters, respectively. 



*K = A , 



(2.6) 



where e° = a/(1 + /3)/(l - /3). 

By causality, the spacetime where two shock waves travel in opposite directions along 
the z coordinate, is described everywhere by superimposing the two shock wave metrics, 
except in the future light cone of the collision. Thus, in the boosted frame the spacetime 
geometry reads 



-dudv + 



vu9{u) ^,,\'^ 



1 



1 + 



vut)(u] 



2-p 



\v9{v) 
2 



1 



dp^ 



2p 



1 



dn^ 



D-3 5 



(2.7) 



which is valid everywhere except in the future light cone oi u = v = 0. As a consequence 
of the collision a strong burst of gravitational radiation (followed by a tail) is expected to 
be emitted. An illustration of this signal in the centre of mass/energy and in the boosted 
frame is presented in Fig. 1. 



3 The perturbative metric in the future of the colhsion 

We shall now compute the geometry in the future of the collision to first order in pertur- 
bation theory. The main idea is that in the boosted frame one shock wave can carry much 
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more energy than the other, since 



V 

A 



.2a 



1 + 

1-/3 



hm — = oo 

p^i A 



Thus, one may face the wave travehng in the —z direction as a small perturbation of the 
wave traveling in the +z direction. Since the geometry of the latter is flat for n > 0, we 
make a perturbative expansion of the Einstein equations around flat space, in order to 
compute the metric in the future of the collision {v^u > 0). To first order this amounts 
to solving the linearised Einstein equations around flat Minkowski spacetime, subject to a 
boundary condition given by the metric (2.7) in the limit n = 0"''. 

The computation will be performed in the more intuitive Brinkmann type coordinates. 
To obtain the boundary condition we transform back to these coordinates on u = O"*". 
Using (2.4) with k replaced by v we find 



dx°' dx^ . 
dxi^ dx' 



:9af5 



dx^ 



dx>^ 



/ 1 



=0+ 



4 
2p 



p 



(3.1) 



where gai3 is the metric for the weak shock by taking u = in (2.7). Observe that the 
dimension of the parameters u and A is [Length] so we choose to perform the following 
rescaling to dimensionless coordinates {u,v,x^) — )■ i'^^^^~^\y/2u, \^v, x^). Performing the 
matrix multiplication and expressing everything in the unbarred dimensionless coordinates, 
we write the geometry, on tt = 0^, as 



llD-Z 



where the flat metric is 



(3.2) 



r]fj,^dx^dx'^ = —2dudv + Sijdx^dx^ . 



(3.3) 



The perturbations^ are 



/^W = {D- 3)^''h{v,p) , = -^V2{D - 3Wh{v,p) , 

/^lj^=2(-5., + (I?-2)^) h{v,p), 



(3.4) 



where 



p 2 



hS^ = ls., + {D-2){D-A)^]h{v,p)\ 



h{v,p) = -— (y2v -<^je (^V2v - . 



(3.5) 



(3.6) 



■^The metric (3.2) is exact on u = 0"^, even thougli it has been written as a perturbation of flat spacetime. 
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Observe that the step function jumps at -y = $/\/2. To understand why, recah that, from 
(2.7), the cohision occurs at u = = in Rosen coordinates. From (2.4), taking into 
account the rescahng to dimensionless coordinates performed above, the coUision takes 
place at 

in Brinkmann coordinates. The future hght cone of the colhsion has two branches: u = 
0,v > and u > 0,v = 0. In terms of the Brinkmann coordinates these conditions read: 

$(p) , $(p) $'(p)2 ^ / u^'(p) 

u = 0, v> — ^ , and u>0, v = — ^ , where p = p 1 + 



In order to understand the interaction, let us follow the null generators of the shock 
with support at D = (which we call the weak shock). Before the collision these can 
be parametrised as 



u = A, v = 0, x' = e ^p= = i . 

To go beyond the collision we look at Brinkmann coordinates. From (2.4): 



n = A , v = e{A) ( + ^ ) , x^ = C\l- ' j . (3.7) 

In terms of these coordinates the focusing effect on the weak shock null generators as they 
cross li = is clear - Fig. 2, 3 and 4. Observe that there is a qualitative difference between 
the behaviour in D = 4 and D > 4: the discontinuous jump in v of the weak shock null 
generators is always positive in D > 4 but it becomes negative for large p in D = 4. This 
will give rise to a different domain in the time integration when computing the overall 
energy emitted. In all cases, nevertheless, the generators will focus, generating a caustic, 
at \/2A = C^^^, corresponding to 



, V2i 



- ln(^/2u) + 1 , D = 4 



D-2 1 



D-A 



D~A 1 



{y/2u)^ 



(3.8) 



After crossing at the caustic the generators enter the curved region of the spacetime. Thus, 
until they cross, the light cone structure is that of the Minkowski spacetime (3.3). It is 
simple to see that the deflection angle a undergone by a null shock generator at distance 
p from the axis obeys, in Brinkmann coordinates, 

tana = — ^— r • (3.9) 

Thus the focusing increases (decreases) with D for short (long) distances, as expected from 
the behaviour of the gravitational force. 
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shock null generators 



Figure 2. Evolution of the weak shock null generators (blue arrows) from the viewpoint of 
Brinkmann eoordinates in the boosted frame in D = A. For m < they are at w = 0; then 
the generators undergo a discontinuity in u at m = 0, which is p dependent and negative for p > 1. 
They jump to the collision surface (green lines, at w = and v ~ — v^lnp), gain shear and focus 
along the caustic (red line). 



Figure 3. Evolution of the weak shock null generators (blue arrows) from the viewpoint of 
Brinkmann coordinates in the boosted frame in Z? > 4. For u < they are at z; = 0; then 
the generators undergo a discontinuity in w at u = 0, which is p dependent but always positive. 
They jump to the collision surface (green lines, at w = 0, v V2/[{D - 4:)p^-'^]), gain shear and 
focus along the caustic (red line). 

3.1 Future development of the metric 

To the future of n = the metric will have a perturbative expansion of the type 





i=l 



(3.10) 




Figure 4. Diagram illustrating (a section of) the spatial trajectories of the null generators of 
the weak shock, exhibiting their focusing after u = 0. Points at the axis (p = 0) will be hit, 
simultaneously, by an infinite number of null generators, two of each are in the represented section. 
After this event, which is at the caustic, spacetime becomes curved therein, and therefore the 
continuation of the rays beyond the axis is merely illustrative and is represented as a dotted line. 
For points outside the axis, such as point V, this diagram suggests an initial radiation signal, 
associated to ray 1, followed by a later burst, associated to ray 2. We shall see this interpretation 
matches the wave forms computed below. 

(i) 

and we should be able to find hiil by successively solving the Einstein equations at each 
order in X/iy, with the boundary condition expressed before. To first order we therefore have 
to solve the linearised D-dimensional Einstein equations around Minkowski space, so in the 
remainder we consider the linear perturbation hfj^i^ in (3.10) and for notational simplicity 
drop the "(1)" label. We denote the trace reversed perturbation by barring it h^^, = 
h^iu — 'nfj,iyh/2, and perform a coordinate transformation of the form — )• rr^^ + {X/v)^^ 

h^u h^u = + '^^{^l,u) ■ (3-11) 
This can be chosen such that the De Bonder gauge condition is obeyed in u > 

/i^"/3_^ = 0, (3.12) 

so that the linearised Einstein equations become simply a set of decoupled wave equations 
in Minkowski space for each component: 

uh^^ = Q ^ ^-2dud, + df)h^, = o . (3.13) 

In [12], and integral solution for the wave equation OF = in n > with initial data on 
u = was used. In appendix A, we show that in D dimensions, the corresponding solution 
is 

F(n, V, X,) = ^ / d^-2x'aj7^F(0, v', x[) , (3.14) 

(27rn) 2 J 
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{u,v,x) 




Figure 5. The intersection of tlie past liglit cone of a spacetime event (u, u, x), to tlie future of 
u ^ = V, witli tlie null hypcrsurfacc w = 0, is a parabola (dashed blue curve); its points (0,w',x') 
obey (3.15). This diagram represents the causal structure of the background metric for unspecified 
initial data. 



where, for each x' , v' defines points at the u = hypersurface, which are on the past light 
cone of the event {u,v,Xi) (cf. Fig. 5): 



V = V 



\x — X 
2u 



/|2 



(3.15) 



and the derivative operator is suitably defined in Fourier space (with respect to v) for odd 
dimensions (see appendix A). Similarly to [12], the gauge condition (3.12) is obeyed in 
u > if 



h 



N 13 



a,l3v\u=0 







(3.16) 



To see this we first contract the wave equation with another derivative to get Oh^u 
It then follows from the integral solution (3.14), that if 



N ,u 



0. 



h 



N IS 



a ,13 







(3.17) 



u=0 



then the De Bonder condition holds everywhere. In particular (3.16) is sufficient, since it 
implies (3.17). 

3.2 Gauge fixing at n = 

Condition (3.16) can be written, in terms of the perturbations in the original coordinate 
system, using (3.11) (on u = 0). Then, using (3.13) and the fact that the trace of the first 
order metric perturbation h vanishes on n = 0, we find, that (3.16) implies, on u = 



1 



h +h " 

'"Ul, V I '"uu, V 



1 



-h- 

2 jv, I 



-{D-3){D-2) 



P 



e [V2v-^j ,(3.18) 
(3.19) 
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where the right hand side of (3.18) and (3.19) was computed using (3.4). We look for a 
solution of (3.18) and (3.19) which has a power series expansion around u = 0: 

W^,v,x^) = eW(^;,x*) +<«(t;,x^) + . . . (3.20) 

One such solution is ^^"^^ = = ^^-^^ = and 

^ ^ J»-3)(I.-2) ^^^^ _ ^p^^^^ ^) ^ (3^^^, 

Applying this gauge transformation, only the uu first order perturbation changes in this 
new (De Bonder) gauge, on u = 0, which becomes 

C = (^-3)$"%,P) + 2d'^ • (3.22) 

Since the metric perturbation is traceless on n = 0, it is traceless everywhere, due to (3.14) 
(because dyh'^ {0,v, x^) = so the solution is identically zero). Thus (3.13) becomes 

nh^^ = ^ {-2dud, + df) h^, = . (3.23) 

The general solution of (3.23) in u > is again obtained using (3.14) with the replacement 
F — )• /i^. Observe that the components of the metric perturbation that vanish at u = 0, 
will vanish when u > 0. We can boost back to the centre of mass frame simply by using 

the replacements {u,v} — ?• {e~°'u,e°'v}. At linear order, we can use xj^ = x^ for the 

1 1 

coordinates in the new gauge. Finally, we can change from units of v^-^ to units of k^-s 
and obtain the linearised, centre of mass frame metric in u > 

-2dudv + h^ui'^^ ^) Xk)du^ + 2/i^(n, v, Xk)dudx^ + \5ij + hfj {u, v, Xk)] dx^dx^ . 



2 

KO-3 

(3.24) 

4 Extracting the gravitational radiation 

To extract the gravitational radiation produced in the collision we shall use the Landau- 
Lifshitz pseudo-tensor [28], which was generalised to higher dimensions in [32]. In terms 
of our perturbations h^j^{u,v,Xi), taking into account that they are traceless, it reads (we 
omit from now on the superscript N for notational simplicity) [32]: 



(4.1) 



Despite not being unique or gauge-invariant it is well known that the integral 

^radiated = / tf^nidSdi , (4.2) 



computed on a 'distant' surface with area element dS outward unit normal n*, is a gauge- 
invariant well defined energy [33]. 
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Following [12], we shall compute the power emitted along the vicinity of the 6 = n 
direction (in the CM frame), corresponding to the negative z direction. This maps to 
the region where perturbation theory should be valid in the boosted frame (cf. Fig. 1). 
So we shall need the flux along the z direction which is given by {t = {v + u) /^/2 and 
z = {v-u)/y/2): 

til = \ (til - til) ■ (4.3) 

The first (second) term corresponds to the flux across a v = constant (n = constant) 
surface. Thus, in order to determine the flux in the 9 = tt direction we need only the 
second term. Then, from (4.1) (observe that the flrst two terms are zero by the De Bonder 
gauge condition): 

The total radiation emitted will be computed, to first order, under the assumption that 
dE/d cos 6 is isotropic (we shall further discuss this approximation below). Thus we in- 
tegrate the power emitted inside the narrow cone around the 9 = tt — 6 = axis. Using 
dS = r^~'^dQD-2, taking the limit close to the axis and multiplying by the area of the 
sphere of radius r, this energy is 

irradiated = . (r^p""-^ [ h'\,hij,„dt) . (4.5) 



32ttGd e-^o, 



As expected, only the ij components (which did not change with our gauge choice) deter- 
mine the energy emitted in gravitational radiation. In appendix B we show that, in fact, 
there is only one independent quantity that determines the integrand in the last equation, 
denoted E = E(u,v,p). In terms of E, we can express the radiated energy as 

-E'radiatcd = ^ hm ( r"^ p^"'^ I {E^yfdt ) = eradiatcdS^U . (4.6) 

dZ'KLjD U - 6 e%0,r^oo \ J J 

We identify, inside the parenthesis, the analogous of Bondi's news function used in [12]. 
The coordinates used in the last formula are dimensionful. The fraction of the energy 
radiated in gravitational waves eradiated defined in (4.6) is given by, 

eradiated = \^ — ^ „ Hm ( / (rp^^,^,)^dt ) , (4.7) 



8 -D — 3 e^o,r-5-cx> 
and, from appendix B 

E, = r ff^^i [41 -4)^] «'^> {< - 4) . (4.8) 

(27ru)^~ -'0 P J-i L ^ 

where 

p^ - 2pp'x + p'^ , _ ^{p') 

v,=v , v, = -^. (4.9) 

All coordinates in (4.7) and (4.8) are now the dimensionless coordinates used in Section 3. 
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Observe that at p = 0, the argument of the delta function has no x dependence and 
the X integral can be immediately performed to yield zero. This can be interpreted as the 
result of a destructive interference phenomenon. To see this consider for instance the /in 
component of the initial data (3.4), which can be written as 



The angular part of the perturbation is a scalar harmonic with £ = 2 (or ^ = — 2 in D = 4) 
on the {D — 3)-sphere. Integrating such scalar harmonic on the {D — 3)-sphere gives zero. 
This is the integral appearing in (4.8) at the axis and may be seen as a cancellation between 
the different phases coming from different points on the integration circle of radius p'. 

The vanishing of the power radiated at the axis in a head on collision of two equal 
objects (like particles or black holes) can also be physically interpreted as follows. Gravi- 
tational radiation emission is determined by the variation of the gravitational quadrupole. 
Taking the observation point at the collision axis, and behind one of the objects, no 
quadrupole variation is observed.'^ 

4.1 Integration limits 

We now discuss the integration domain for the time integration in (4.7). Equation (4.8) 
suggests that the gravitational radiation seen at a spacetime point V to the future of the 
collision, originates from the points on the collision surface wherein the delta function has 
support; these are the points at the intersection of the past light cone of V with the collision 
surface, as expected.^ 

Instead of using the spacetime coordinates (u, v, p) it is convenient to introduce a 
retarded time coordinate r = t — r, where the radial coordinate is r = \J ^ p^ and an 
angular coordinate ^ by p = rsin^. Thus we wish to know the first instant t\ for which 
the point specified by (r, Q) receives gravitational radiation. This event occurs when the 
intersection of its past light cone with u = becomes tangent to the collision surface and 
is determined in a similar way to the determination of the caustic in Section 3. One finds 
that Ti = ri(r, Q\ together with the auxiliary variable p, is obtained by solving the system 
of equations 



for specified p and d and s = +1. The intersection of the past light cone of points (ri, r, 0) 
(points 1 to 3 in Fig. 6) with m = is plotted in Fig. 7, for = tt and Fig. 8, outside the 
axis. As claimed these curves are tangent to the collision surface. 

^We thank U. Sperhake for this observation. 

^Observe, however, that an integration in v was aheady performed to get (4.8). For even D this 
integration is performed using a delta function which enforces the support of the integrand on the light 
cone. But for odd D the initial integral has support also inside the light cone, cf. Appendix A. 



f 2x2 - 1 



L» = 4 




(4.10) 




(4.11) 
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Figure 6. The caustic is represented in a spacetime diagram for the case D > A. Along the caustic 
three events 1 — 3 are selected. For event 3 its past light cone is drawn and its intersection with 
the 7i = surface is represented as the blue dashed line (parabola). This is tangent to the collision 
surface u = and v = $/\/2, represented by the solid (green) lines. The generators of the shock 
traveling along u that emerge from the intersection points will focus and converge at 3. For an 
observer at fixed z, following the worldline represented by the r axis, no gravitational radiation 
will be observed before r = ti. In Fig. 7 (right panel) the collision lines and the intersections with 
u = of the past light cones of other points along the caustic are represented. 




Figure 7. Collision surface (at it = and v = $/-\/2 - green solid lines) and the past light cones 
at w = (blue dashed lines) of various spacetime points along the caustic for D = A (left panel) 
and D ^ 5 (right panel). The past light cones are tangent to the collision surface. The latter case 
is qualitatively similar to any D > 5. The numbering 1-3 in the right panel corresponds to the 
spacetime points with similar numbering in Fig. 6. 

Thus the initial time for observing gravitational radiation, ti, is positive, for D > 4 
and goes to zero as r goes to infinity (cf. Fig. 3). By contrast, in D = 4, ri becomes 
negative and approaches negative infinity as r goes to infinity (cf. Fig. 2). Thus, the time 
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Figure 8. Same as the previous figure but now for spacctime points away from tire axis (witlr tlie 
same u = 1 but different distanees from the axis p = 0, 1, 3) such that their past light eone at m = 
is tangent to the collisions surface, for _D = 4 (left panel) and _D = 5 (right panel). The earliest 
intersection always occurs for x = 1. 

integration will have different domains in = 4 and D > 4. 

The previous computation of ti amounts to computing the retarded time at which 
ray 1 intersects point V in Fig. 4. The retarded time, T2, for which ray 2 intersects V is 
computed by solving (4.11) with s = — 1. This coincides with a second peak in the wave 
forms exhibited in the next section. 

4.2 Numerical evaluation 

In this section we finally obtain the wave forms for several dimensions and integrate the 
radiated power to estimate the amount of gravitational radiation emitted in the collision. 
First we perform the angular integral in (4.8) (see appendix B where the polynomials of 
degree M + 2, p(*^+2) and are defined), to obtain 

M+i^ r f dp' p(A/+2)(^^ 



in odd dimensions and 



in even dimensions. T) is defined such that 

2pp' 

where U = t + 2rsin^(0/2), T = t + 2r cos'^{6 /2) — p^ /U and we are now expressing the 
result in coordinates {r, r, 0} as in Eq. (4.11). 

In the remaining numerical analysis, we present only results for even D. Concerning 
odd D, we have attempted to evaluate (4.12); the wave forms obtained (numerically). 
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r — ri T — Ti T — Ti 



Figure 9. Wave forms: These plots show the even D wave forms (4.13) for fixed 9, close to the 
axis of symmetry and for various r. Note the initial step which is due to the impulsive nature of 
the initial conditions associated with the Aichclburg-Scxl shocks. This step is exactly at t — ti, 
whereas the sharp peak is at t = T2, for each curve. 

however, contained non-integrable tails when squared to get the radiated power. Our best 
hint, at the moment, is that these divergences are similar to those present in the second 
order computation [13]. Nevertheless, our aim is to investigate the variation of the radiated 
energy with D. Thus the even case already shows the general trend. It seems plausible 
that the result for D odd will interpolate between the even D cases we shall exhibit. 

To evaluate (4.13), we wrote firstly a test code in mathematica7 and then a code in 
the CH — h language, using the numerical integration routines of the Gnu Standard Library 
(GSL). The purpose was to check the two codes in some cases and use the (lower level) CH — h 
code for generating all the data in a practical amount of time. The domain of integration 
P was determined by looking at the roots of two polynomials constructed from and 
using the root finding routines in MATHEMATICA or the GSL library. In general, depending 
on the spacetime point, we found either an empty integration domain (giving a value of 
exactly zero to the wave form), or a non-empty domain which can be simply connected or 
the union of two disconnected domains. The integrable singularities at the points x^, = ±1 
where removed explicitly through a change of variable around each singularity. In general, 
for the wave forms, we have demanded a relative error of 10~^, and 10~^ for the final r 
integral of the radiated power. 

In Fig. 9, we present some wave forms {D even) as a function of r. Since we are 
interested in the limit near the negative z axis (where our approximations are justified), 
we use examples with 6 away from the axis by 0.01 radians. A first observation is that we 
have shifted r by ri for all wave forms such that zero corresponds to the beginning of the 
radiation burst. A consistency check is that, indeed, the numerically determined integration 
domain agrees exactly (within very small numerical errors) with this expectation. Thus 
all the wave forms have a sudden step at zero which arises naturally from the numerical 
code. Another feature which is verified numerically, is that the (integrable) singular peak 
of radiation observed in all plots, occurs exactly at r = r2. This is explicitly exhibited in 
Fig. 10, where besides the shift of r by ri we have rescaled the axis by Ar = T2 — ti. 
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Figure 10. Wave forms in a rescaled r coordinate: For all wave forms, the signal jumps to a 
non-zero value at r = ri and peaks at r = T2. These are the values of r at which rays 1 and 2 of 
Fig. 4 reach the observation point, respectively. 



An advantage of our numerical method compared to [12] is that no approximations 
were made in the integral. In fact in [12], the limit r — )• +00 was taken inside the integral 
before actually performing the integration. This removes the p'^ term in (4.14), which 
is responsible for excluding regions of the domain of integration with large p'. In four 
dimensions, due to the logarithmically growing <I> profile, it turns out that such regions of 
large p' are indeed cut off from the integral by the In p' term. Thus in the large r limit, the 
wave form does approach the one provided in Fig. 4 of [12]. This can be seen in the first 
plot as r becomes large. However the approximation in [12] makes it less clear when the 
radiation burst starts and peaks. In our approach we simply use the full numerical integral 
without further approximations, for finite r, and find the limit numerically. Finally, for 
D > 4, the procedure in [12] is not valid, because it would amount to neglecting the only 
term growing with p' in the definition of x^, so there would be no regularisation of the 
domain and the integral would not converge. For D > 4, there is a notable difference as we 
take the large r limit, namely that the radiation pulse becomes more concentrated around 
r = T2. In fact, the same happens in four dimensions, however more slowly (as r increases) 
due to the logarithmic nature of the profile of the shock wave ^. The other difference is 
that the wave form acquires one more oscillation in each jump to the next D even. Finally 
we note (without plotting explicitly) that the effect of decreasing 6 (with r fixed), is the 
same as when we decrease r with 6 fixed®, that is, the radiation burst becomes sharper and 
occurs in a smaller period ~ At as we move closer to the axis. 

After we obtain the wave forms for a certain r,6, we can compute the estimate for 
the fraction of radiated energy, by computing the integral in (4.7) and taking the limit 
r — )• +00, 9^0. In Fig. 11 we show plots of the radiated fraction of energy as a function 
of r for several 9. We have performed a fit of each numerical curve to extract the limit, by 

^This is clear by looking at the ranges of r for the different curves in each plot where, in D = 4, r varies 
by several orders of magnitude, whereas for D > 4 it only varies roughly by one order of magnitude and 
the separation between curves is similar. 

®This should not be surprising since the r and 9 dependence appears mostly in p = rsin^ factors. 
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Figure 11. Limiting fractions: These plots show e(r), whose Hmit for r large and small n — 9 
gives the estimate Eradiated- The best estimate for the limit can be read from the constant term in 
the asymptotic form used to fit the numerical points of the red curves (for which the angle is the 
smallest). 

using an asymptotic expansion of the form e(r) ~ eradiated + a/r (indicated next to each 
curve). Since we have computed this function for very large values of r we obtain a fit 
which is basically indistinguishable from the numerical tails, so eradiated is extracted with 
a relative error of less than 0.001. Note that the difference in the shape of the curves for 
different 9, is again smaller for the D = 4 case due to the logarithmic dependence of the 
shock wave profile. Since we have used very small angles, all the asymptotic values in the 
different curves are consistent with the limiting values quoted in the introduction (within 
the relative precision of 0.001 and the error from the a/r correction). 

4.3 Discussion of the method 

The computation performed in this paper relies on perturbation theory. Why should 
perturbation theory be valid in describing a non-linear effect (like black hole formation) and 
an associated strong field effect, such as the generation of a strong burst of gravitational 
radiation? In [12] this is justified as follows. Mathematically, the boost provides a small 
parameter for expanding the geometry (the ratio of energies of weak to strong shock waves). 
Analysing (3.4), one concludes that in the boosted frame the initial conditions are indeed 
perturbative except in a small vicinity of /? = and sufficiently large v > 0. This has a 
physical picture. The null generators of the weak shock wave (in the boosted frame) are 
not only bent towards the axis of the collision, as they cross the strong shock, but they 
also suffer a redshift which becomes larger as the null generator approaches the axis and 
infinitely large for generators at the axis. This is manifest in Fig. 2 and 3. The radiation 
that hits a spacetime point (not at the axis), therefore, comes first from the far field region 
of the collision (ray 1 in Fig. 4), which is given a large head start as compared to its near 
field counterpart (ray 2 in Fig. 4). Since gravity is weak in the far field region, perturbation 
theory should therefore be accurate in describing this part of the signal. Moreover, due to 
focusing, the overall amount of radiation measured near the axis, at very large distances, 
coming from the far field region need not be small. Therefore, it is concluded in [12] that 
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this problem provides one example where perturbation theory can tackle a strong field 
effect. 

The interpretation of the wave forms we have provided in this paper shows, however, 
that the above argument must be considered with care. The wave forms exhibited in Fig. 
9 and 10 are dominated by the signal near T2, which is the signal that crosses the near 
field region. This indicates therefore, that there is a large contribution to the signal from 
outside the region where linear theory is clearly applicable and a correction from higher 
order terms in perturbation theory is expected. But the near field region experienced by ray 
2 becomes smaller at very large r and small 9. Furthermore, as we have learned in D = 4, 
the first order perturbation provides a better estimate than the trapped surface argument 
and, more importantly, the second order perturbation provides a result consistent with the 
numerical simulations. We find this an important motivation to pursue this computation. 
It will be interesting to see if the agreement between numerical results and perturbation 
theory still holds in higher dimensions. 

In order to compute the total radiation emitted, we have determined the radiation 
in the vicinity of = vr and extrapolated it under an isotropy assumption for dE/d cos 6, 
following [12]. This is partly justified by the Zero Frequency Limit approximation [34], 
which predicts an isotropic angular distribution of gravitational radiation (with two blind 
spots along the symmetry axis) for the collision of two equal mass black holes as the speed 
of light is approached. Nevertheless, the second order computation of D'Eath and Payne 
introduces an angular dependence [14], which seems to be crucial for a more accurate 
estimate. 

One final remark concerning odd D. It is well known that the Green's functions of 
the D'Alembertian operator have support on the light cone in even D and both on and 
inside the light cone for odd D, cf. (A. 12) and (A. 13) (see, e.g. [35-37]). This leads to 
the unfamiliar property that odd D dimensional Minkowski space behaves as a dispersive 
medium. It may seem that the integral solution (3.14) is missing this property for odd D, 
since it has support solely at the intersection of the past light cone with the initial data 
surface (as illustrated in Fig. 5). But it is not so. The equivalence between the integral 
solution and the Green's function method is shown in appendix A. 

5 Conclusions 

In this paper we have provided an estimate for the energy radiated into gravitational waves 
in a collision of shock waves in a D dimensional spacetime. These shock waves provide an 
intuitive description, due to the infinite Lorentz boost, of the gravitational field of ultra- 
relativistic particles. Moreover, they do capture their gravitational interactions. Indeed, 
the scattering amplitude for two (scalar) particles in the eikonal regime of perturbative 
quantum gravity [38], can be equivalently computed by analysing the scattering of a (test) 
plane wave in a shock wave background [11] (see also [39]). Thus our computation applies 
to particle collisions in a regime wherein their interaction is dominated by gravity, and well 
described by classical general relativity. This is the transplanckian scattering regime [11]. 
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If the fundamental Planck scale is of the order of the TeV scale, a key issue for phe- 
nomenology is how good the classical and semi-classical approximations are to describe the 
black holes that could be produced in particle collisions. In other words how much of the 
initial centre of mass energy of the process stays in the final black hole. The answer essen- 
tially amounts to understand how much energy is lost in gravitational radiation. The best 
estimates for this radiated energy, to date, come from apparent horizon computations [16]. 
In four dimensions the apparent horizon estimates are off by a factor of two as compared to 
improved values obtained either from perturbative computations in colliding shock waves 
backgrounds [12-14] or numerical relativity results [21]. The final black hole is therefore 
more massive (and hence more classical) than what has been anticipated by apparent hori- 
zon estimates. A similar conclusion in higher dimensions would have phenomenological 
relevance. 

In this paper we have computed the metric in the future of the collision of two shock 
waves in first order perturbation theory and extracted the corresponding gravitational 
radiation. Although our result should only be faced as an estimate, what we know from 
D = 4 indicates that it is an improved estimate as compared to that obtained from the 
analysis of trapped surfaces [16]. 

The apparent horizon estimate gives the same trend for the D dependence as our 
result. This trend is also observed in numerical relativity results for a low energy collision 
of black holes in D = 4 and D = 5: these show that the energy radiated increases, from 
D = 4 to D = 5 [30] . Curiously, another bound that applies to these low energy collisions 
exhibits an opposite trend. From the area law, that bound may be obtained, as first shown 
by Hawking [40] , for the energy emitted into gravitational radiation in a collision of non- 
spinning, equal mass black holes starting from rest at infinite distance. This bound starts 
at precisely the same value, 29.3%, as the one discussed herein for high energy collisions, 
but it decreases with dimension. Another context where a similar decrease with D occurs, 
though in a different regime - that of extreme mass ratio, was studied in [41-43], where the 
gravitational radiation emitted by point particles falling into a black hole was computed. 

This paper also clarifies two aspects of the first order computation by D'Eath and 
Payne [12]. The first is the approximation procedure in determining the emitted power, 
namely neglecting p'^ terms in the computation of (4.8). Since we have evaluated (4.8) 
numerically, no such approximation was used, and our result is consistent with that of 
D'Eath and Payne within less than one percent. This shows that the approximation used 
in [12] is robust in D = 4. The second is the time integration domain. It seems strange, 
at first sight, that there is a non-vanishing gravitational signal for all r G M in D = 4, 
in the limit r — t- oo for the observation point. This is a consequence of the nature of the 
coordinate transformation from Brinkmann to Rosen in D = 4 and the jump of the null 
generators of the weak shock to negative values of v as they cross the strong shock (cf. 
Fig. 2), for large p. In D > 4 the coordinate transformation is different and so is the 
time integration, which has support only for r G M"*". Moreover, the observation we have 
made concerning the matching of the peaks in the obtained wave forms with a simple ray 
analysis, is useful in understanding the shape of the wave forms, which was not discussed 
in [12]. 
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Finally, we hope that this paper will help setting the stage for further use of this tech- 
nique, since it can be applied to various generalisations of the problem we are considering. 
In particular, one immediate goal is to consider the second order in perturbation theory. 
We observe that the initial conditions are exact to second order which might be at the basis 
of the good agreement between the perturbative computation and the numerical relativity 
result in D = 4. 
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A General integral solution 

In this section we generalise the integral solution found in [12] to higher dimensions. The 
most standard method is to use Green's functions. In our case, however, there is a simple 
way of finding the solution by using an ansatz inspired by the four dimensional case. First 
we provide such a derivation and after we check it using the well know Green's functions 
method in Z)-dimensional Minkowski spacetime. 

A.l Integral operator ansatz 

Before attempting to generalise the integral solution of the wave equation OF = in n > 
with initial data at u = provided in [12], let us start by investigating how it works in 
four dimensions. In that case the solution is 

F{u, V, = ^ / ^'^'^^(0' ' 
where, for each x' , v' defines points, at u = 0, on the past light cone of the event (u, v, Xi): 

y' = v- . (A.2) 

2u ^ ' 

We can compute partial derivatives of F by acting directly on the explicit u dependence, 
and the implicit dependence on u,v,Xi through v' . Then using 

d d d dv' d \x — x'\'^ d d dv' d x* — x'* d 

dv dv' ' du du dv' dv' ' dxi dxi dv' u dv' ' 

(A.3) 

we can compute OF and check it is indeed zero. Furthermore, we need to check the initial 
condition. If we try to take the limit n — t- directly in (A.l) we realise that there is no easy 
way to deal with the v' derivative. To remove it we can work in Fourier space by taking 



- 21 - 



the Fourier transform with respect to v (we denote it F and replace the corresponding 
argument by q) 



I / 12 



F(u, q, X,) = J d^x'^—^^^F{0, q, x\) = j d^x'6,{xi - x'^F{0, q, x'^) , (A.4) 

where we have identified the parameter e = —iu/q and a function 5^ which becomes the 
2-dimensional Dirac delta distribution in the e — t- limit. So, in Fourier space, the initial 
condition is obeyed and the same holds after Fourier inversion. 
In higher dimensions we try a similar form 

Fiu, v,Xi) = j^l d^-2x'O,.F(0, v', x'i) , (A.5) 

where for the moment we leave the overall normalisation N arbitrary, v' is as before, and 
a is to be determined as well as the O^' operator acting on F. We can again compute 
□F = to obtain 

(D-2- 2a)^^^ j d^-^x'd,,O,,F{0,v',x';) = . (A.6) 

Thus this integral ansatz must have a = {D — 2) /2. Incidentally, this is exactly the correct 
number of u powers needed to obtain a, D — 2 dimensional delta function by repeating the 
reasoning for the initial condition. Again, we Fourier transform with respect to v to obtain 



- —iq- 



F{u,q,Xi) = -^^4?^ / d^-^x'- / dv'e'i'''O^:F{0,v',x[) 



N 
(2vr) — 



2 



So in the n — t- limit we obtain 



(2?™) 

d^~^x'S,{x-x'){-iq)-^ I dv'e''"''O^>F{0,v',x'i) . (A.7) 



D-2 

F{0,q,Xi) = ^^4?^ ^jrr / dv'e^'^^' O,,F{0,v' ,Xi) . (A. 



N 



-iqY 



2 



In higher dimensions, for the right hand side to reduce to the left hand side, we need a 
generalisation of the Fourier transform of the operator 0^' — )• d^, . This is defined 
through its Fourier transform when acting on a function f{v) (so it is a distributional 
operator) 

'e''^^0,f{v) = {-iq)'^f{q). (A.9) 



I' 



Using (A.9) and choosing N = (27r)^^~^^/^, then (A. 8) is solved. In even dimensions Oy is 
just a partial derivative. However in odd dimensions we have a fractional partial derivative, 
thus formally, the general integral solution is 

F{u, V, Xi) = [ d^'-^x'd^FiO, v', x'i) . (A.IO) 

(2™) 2 J 
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A. 2 Green's function solution 



A more standard derivation, consists of reducing our Cauchy problem with initial conditions 
on the null hypersurface w = 0, to a wave equation problem for a distribution with a 
source term. Then the general solution in D dimensions can be found, for example, from 
Theorem 6.3.1 of [44] 



F{u,v,Xi) = 2 
The Green's function for even D is 

G{u,v,Xi) = 



du'dv'd^^^x' G{u -u',v- v', x^ - x'^5{u)d^,>F{0, v',x'^) . 



2{2ttu) — 



\x\ 
'2u 



(A.ll) 



(A.12) 



where the superscript denotes denotes the distributional derivative of order {D — 2)/2, of 
the Dirac delta distribution. For odd D 



G{u, V, Xi 



2u 



D-2 
2 



(A.13) 



where this distribution is to be understood as giving the finite part of the integral when 
acting on functions. In even dimensions, we can perform the u', v' integration immediately 
and use the action of the distributional derivative of the delta function to get exactly (A. 10). 
For odd D, we first note that by working in Fourier space, we can define the negative order 
of the delta distribution, (p < 0) as 

1 



5(P)(a 



r(-p)- 



-p-i 



e{x). 



(A.14) 



This definition can be used to extend (recursively) to positive orders, by acting with deriva- 
tives on 6^P\x). For example, the fractional derivative of the delta function of order 1/2, 
is simply 



d_ 

dx 



5^-^'^\x) 



d 

dx 



\X\ 



TT X 



e{x) 

2x 



5{x) 



.(A.15) 



and any higher order fractional derivative can be defined similarly. We note that up to 
singular delta function terms which we must discard (due to the finite part prescription for 
its integral with a function), we can write the D odd Green's function exactly as (A.12) with 
the fractional order delta distribution. Thus the same integral solution, with a fractional 
partial derivative, holds in odd dimensions in agreement with the first derivation. 



B Solution for the transverse components 

In Fourier space (with respect to f), the solution of (3.23) for the spatial components of 
the metric perturbations reads 

mu,q,Xk) = ^ / d^-2xV''^(-i(?)^/ii}(0,g,x'J . (B.l) 

(2™)— J 
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For convenience we write the initial perturbation as to separate the v dependent term: 

hf^{<d,v,xu) = [5.,, -{D- 2)^) ^ (V2v - e {V2v - cd) . (B.2) 
Then the Fourier transform of the initial perturbation is obtained as 

(27ru)— J V P J P 

(B.3) 

The integration measure and p' are invariant under rotations. So we define a rotation, such 
that Hi = RijXj, is aligned with the yi direction and has magnitude p. Similarly y'^ = Rijx'j, 
so then 



hfAu,q,Xk) 



.-i.-iV2H^r / ,M.\^\p')A ''~'-'r'''' ^^-^ 

R.^ ^.rn^^-o^ J d y ^<^.„ - {D - 2)-^J —^e 

(B.4) 

where 9' projects the vector y' onto the y[ direction. In D = 4, 6' £ [0, 2tt] whereas for 
-D > 4 6*' S [0, tt]. The parity of the integrand with respect to 7/21 •• • ■•y'D-2 always well 
defined (there is only no definite parity with respect to y'^ due to the exponential). Thus 
only the £ = m components are non-zero, either due to the 6im or because if i ^ m, then 
y'^y'm is odd when integrated over y2, ■ ■ ■ ,y'£)^2- Changing to hyperspherical coordinates 

mu,q,x,) = R-^RJ^X3Z1S^ (B.5) 

[2ttu) 2 

/ dp'^'ip')p'^-^ de'isme')''-^e'\ ^" J dno-A 



The angular integrals are obtained straightforwardly using standard integrals on the [D 
2)— sphere, so we obtain 



g-6 

2 



hij{u,q,Xk) — R^^ Rj^ n-2 

{2ttu) 2 



/riT j / p^-2pp'cosfl' + p'^ I f-jp') I 

dp'^'{p')p'^-^ j de'{sm9' f-^e \ ^ /17d-4G^™(cosO , (B.6) 



with 



G,„(cos e') = b,m - diag | (I? - 2) cos^ Q\ ^--| sin^ 0', . . . , sin^ 0' j , (B.7) 

where we have defined J7o = 2 to account for the fact that Q' is a polar angle in D = 4. 
Note that the first angular integral (corresponding to y'^y'i) is just an area factor and the 
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other components can be obtained using the vanishing of the trace. This last property 
is satisfied by the final result (B.7) which works as an independent check of the angular 
integrals. Fourier inversion then yields 



hfj{u,V,Xk) = R,/Rjm 



1 



D~2 
I 2 



{2-KU 



(B.8) 



2u y/2 



The exponent of the delta denotes the order of the derivative, which can be fractional 
or negative as in (A. 15). Note that due to the vanishing of the trace, actually only one 
component of the ij is independent. If we factor out the rotations, we can define the 
rotated 11 perturbation as 



hu{u,v,p) = E{u,v,p) 



(27ru) 



D-2 
2 



J dp' p"'-'t'(/)x 



Jo \ 2ii ■y/2 / 

(B.9) 

Now if we define x = cosO', take a v derivative and use the scaling properties of the 
derivative of the delta distribution 



with 



D-4 



{27rp) — {D-l)Jo p'- 



dp 



(l-x^)Vj , (B.IO) 



U<^> (p') + p'^ - UT 



(B.ll) 



where U = \/2u and T = \f2v — jU . To perform the integral, we can always integrate 
by parts M = [{D — 4)/2] times (the brackets indicate the integer part) without obtaining 
any boundary terms^, so we end up with 



VsnD-ii-iV'^ /■+°^ dp' ^ d^'^ 



1 



(27rp)^~(£» - 1) JO iJ': 

where g = for D even and <? = 1 for D odd. Then in even dimensions (note T) is defined 
such that — 1 < < 1) 

^/8^^z)-4(-l)'^ 



Em 



dp' r 



(27r/))A^+i(D - 1) p>M+2 
and in odd dimensions, using (A. 15), 



fl-x2)^^+l 



(B.13) 



E_, 



D-4l 



■1) 



M+1 



{27rp)^^+i^{D -I) Jo p 



dp' 



1 



[(1 



2\M+2-\ 



dx 



' (2; - X*) 



(B.14) 



^This is because terms of the form (1 — x^)'' are always present (evaluating to zero at the boundaries) 
and the delta functions only have support on the boundary for a subset of points with measure zero. 
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The X integral can be performed by expanding the polynomials and integrating by parts. 
Then if we define the following polynomials 



M+2 



[(1 - ^r'^'] - E 




(B.15) 



fc=0 




(1 _ ^2)M+|J ^ ^ d^^k 
k=0 



(B.16) 



^^M+2 [ 




(B.17) 



we obtain (4.12) in odd dimensions and (4.13) in even dimensions. 
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